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Method For Processing Raw Magnetic Resonance Data, 
Magn tic Resonance Tomography Unit And Magnetic 
R sonance Spectroscopy Unit 

[0001] The present application hereby claims priority 
under 35 U.S.C. §119 on German patent application 
number DE 10247150.9 filed October 9, 2002, the entire 
contents of which are hereby incorporated herein by 
reference . 

Field of the Invention 

[0002] The invention generally relates to a method for 
postprocessing raw magnetic resonance data. The 
invention also generally relates to a magnetic 
resonance tomography unit or a magnetic resonance 
spectroscopy unit which is respectively matched to such 
a method for postprocessing raw magnetic resonance 
data . 

Background of the Invention 

[0003] The aim of processing raw magnetic resonance data, 
which can originate both from magnetic resonance 
spectroscopy units and from magnetic resonance 
tomography units, is to extract the medically relevant 
information from the raw magnetic resonance data as 
well as possible. In the case of magnetic resonance 
tomography, the aim is accordingly a picture with high 
resolution of detail and low noise. That is to say, 
the aim is for a picture with high edge sharpness and a 
high signal to noise ratio (SNR) . 

[0004] In the case of magnetic resonance spectroscopy, 
the aim is to obtain from the raw data a linear 
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spectrum with high resolution and with a high SNR. The 
quality in the processing of raw magnetic resonance 
data is thus firstly in the edge sharpness and image 
sharpness, i.e. in sharp pictures with high contrast. 
Secondly, it is in the production of a high SNR, so 
that the essence of the picture or of the spectrum is 
highlighted. 

[0005] In the case of the various approaches to achieving 
this aim, it is usually not possible to increase the 
image sharpness and the SNR simultaneously. This 
drawback has a direct effect on the quality of the 
magnetic resonance tomography and spectroscopy. To 
improve the resolution of detail, radio- frequency 
"blue" noise is superimposed on a picture, for example. 
This produces a subjective increase in the resolution 
of detail. A likewise subjective improvement in the 
resolution is brought about in a similar manner by the 
application of high pass filters to the image data. 
This also amplifies the image noise, however. 

[0006] To improve the SNR in the case of magnetic 
resonance pictures, known filters (such as the Hanning 
filter, the Fermi filter or the cosine filter in the 
frequency domain) are applied to the raw magnetic 
resonance data. The action of such filters is known 
from the literature, e.g.: F.J. Harris, Proc . IEEE, 
Vol. 66, No 1 (1978): "On the Use of Windows for 
Harmonic Analysis with the Discrete Fourier Transform" . 
It is also possible to apply postprocessing filters in 
the space domain, for example a mean filter, median 
filter or ARMA filter. 
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[0007] in addition, US patent 4,463.375 discloses a 
method for medical image processing. This method 
specifies a reduced-noise version of a first processed 
image which has been obtained from a multiplicity of 
measurements in a multipicture system, e.g. using a 
computer tomography unit. The method first involves 
producing a second image with a high SNR from the 
multiplicity of measurements. The first image is then 
processed using a filter which reduces the noise. The 
second image is processed using a filter which is 
complementary to this filter. The weighted combination 
of the two images produces the reduced-noise version of 
the first processed image. 

SUMMARY OF THE INVENTION 

[0008] An embodiment of the invention includes an object 
of specifying a method for processing raw magnetic 
resonance data and a magnetic resonance tomography unit 
and a magnetic resonance spectroscopy unit which 
produce pictures or spectra with high edge sharpness 
without impairing the image noise, i.e. the SNR, in the 
process . 

[0009] For the method, an embodiment of the invention 
achieves an object by virtue of a method having the 
method features below being carried out in order to 
postprocess raw magnetic resonance data. In a first 
step, the raw magnetic resonance data are filtered 
using a first filter. These filtered data are Fourier 
transformed in a second step. The absolute value for 
the Fourier transformed data is then formed in a third 
step and in this way a first magnetic resonance signal 
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is obtained. In a fourth step, the unfiltered, original 
raw magnetic resonance data are Fourier transformed. 
This is followed in the fifth step by second absolute 
value formation for the unfiltered Fourier transformed 
data. Finally, weighted combination of the two magnetic 
resonance signals is performed and in this way the 
postprocessed magnetic resonance signal is obtained. 

[0010] The method includes the weighted combination of 
two magnetic resonance signals, the magnetic resonance 
signals having been obtained from the same raw magnetic 
resonance data once with and once without filtering. 
Through suitable choice of the filter and of the 
parameters for the weighted combination, a picture 
which is optimized in terms of the representation 
wanted from a medical point of view is extracted. 

[0011] One advantage of the method is that the SNR is 
significantly improved following the weighted 
combination of the two magnetic resonance signals and 
that an edge - emphasized magnetic resonance signal is 
produced at the same time. In comparison with the 
method from US patent 4,463,375, the involvement for 
this method is significantly reduced. This occurs 
since it is firstly based only on the raw data from a 
single measurement, and since it is secondly 
independent of an image with a high SNR which is 
produced from a plurality of measurements. 

[0012] In one particularly simple extension of the 
method, the raw magnetic resonance data are filtered 
using a second filter before the second magnetic 
resonance signal is produced. This allows the second 
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magnetic resonance signal to be influenced and thus 
allows it to improve the weighted combination. The use 
of the unfiltered raw magnetic resonance data in the 
original method corresponds to the application of the 
identity filter and is faster to execute. 

[0013] A weighted combination within the context of an 
embodiment of the invention is in this case a 
mathematical combination of the absolute values of the 
two magnetic resonance signals . These have the same 
local coding and represent the echo signal from a 
magnetic resonance measurement, for example in a 
one-dimensional or multidimensional pixel structure, 
following the Fourier transformation. In this case, the 
combination, which is usually addition, is spatially 
consistent, i.e. the same spatial areas are combined . 

[0014] By way of example, combination takes place pixel 
by pixel. The weighting is effected using a weighting 
factor. For a spatial area of one of the two magnetic 
resonance signals, this weighting factor depends on the 
information content, e.g. the intensity, of the same 
spatial area in the other magnetic resonance signal . By 
way of example, the weighting factor can have a linear 
dependency on the information content or can have any 
other beneficial nonlinear dependency on the 
information content of the spatial area. The weighting 
factor can be produced, by way of example, by a 
mathematical step function which sets the weighting 
factor to a small value or zero below a previously 
stipulated threshold value for the absolute value of 
the magnetic resonance signal in question. 
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[0015] One advantage of weighted combination is that the 
SNR is significantly improved following the weighted 
combination of the two filters. Further, this data 
weighting produces an edge - emphasized profile in the 
magnetic resonance signal. 

[0016] The method relates to the postprocessing of raw 
magnetic resonance data from a single measurement and 
is not dependent on whether the raw data belong to a 
single picture or whether they belong to a picture from 
a series of pictures which show a development over 
time, for example. The raw magnetic resonance data can 
be raw data from a measurement in a magnetic resonance 
spectroscopy unit or else raw data from a picture in a 
magnetic resonance tomography unit. 

[0017] Accordingly, an embodiment of the invention 
achieves an object for a magnetic resonance tomography 
unit by use of a magnetic resonance tomography unit 
which is matched to a method for postprocessing raw 
magnetic resonance data with the method features above. 
In this case, the raw magnetic resonance data are 
referred to as a spin echo signal or an echo signal, 
for example. The magnetic resonance signals correspond 
to the one-dimensional or multidimensional magnetic 
resonance pictures, for example. 

[0018] An embodiment of the invention achieves an object 
for a magnetic resonance spectroscopy unit by use of a 
magnetic resonance spectroscopy unit which is matched 
to a method for postprocessing raw magnetic resonance 
data with the method features above. In this case. 
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spectroscopic data, e.g. the FID (Free Induction Decay) 
signal, are postprocessed. 

[0019] In one particularly beneficial embodiment of the 
method, the filters are applied to the raw magnetic 
resonance data after demodulation thereof. This allows 
the filters to be applied in their customary form of 
representation in the frequency domain. Applying the 
filters after demodulation has the advantage that the 
filters can be very narrowband. 

[0020] In one particularly advantageous embodiment, the 
first filter is a low pass filter which filters rapid 
changes out of the raw data and thus suppresses radio- 
frequency oscillations. By way of example, the low pass 
filter can be a type of Manning filter. 

[0021] In another embodiment the second filter is a high 
pass filter which passes only the rapid changes in the 
signal, which are also caused by noise, for example. 
Magnetic resonance signals which have been high pass 
filtered show an excessive increase in rapidly rising 



[0022] In one particularly advantageous embodiment which 
uses self-weighting for the combination, the weighted 
combination is effected as follows: 



edges . 




[Equation l] 
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[0023] In this case, A and B are magnetic resonance 
signals which are subjected to weighted combination by 
low pass filtering (A) and high pass filtering (B) , for 
example. Amax is the maximxim of the magnetic resonance 
signal A. The two parameters X and k determine the 
contribution level of the magnetic resonance signal B 
to the postprocessed magnetic resonance signal C. The 
quotient of B and A^ax brings about amplified correction 
of the pixels which have a high absolute value in the 
magnetic resonance signal B. The values of X and K need 
to be ascertained empirically on the basis of their 
influence and their optimum effect on the magnetic 
resonance signal. In this case, the values of X and K 
are preferably between 1 and 3 . 

[0024] One advantage of the quotient which normalizes B 
to Ainax is that an uncontrolled signal contribution for 
B is avoided in regions in which there is actually 
nothing to correct or to improve. 

[0025] The weighted combination of suitably filtered 
magnetic resonance signals has the advantage that, in a 
single magnetic resonance measurement, the noise in the 
resultant magnetic resonance signal is lower than with 
linear addition or subtraction of the two filtered 
magnetic resonance signals. 

[0026] In particular embodiments, the method from the 
invention can be applied to raw magnetic resonance data 
which cover a one-dimensional or multidimensional space 
which is to be examined. Depending on the application, 
one -dimensional sections , two-dimensional sectional 
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views or three-dimensional volume representations are 
shown, for example, in magnetic resonance technology. 
The raw magnetic resonance data contain the dimension 
of the respective picture. 

[0027] In one advantageous embodiment, the dimensionality 
of the Fourier transformation performed during the 
method is chosen such that the dimension of the raw 
magnetic resonance data is retained. Thus, the method 
can be applied, by way of example, to two-dimensional 
magnetic resonance sectional views or to 
three-dimensional volume representations. Even in the 
case of magnetic resonance spectra, one-dimensional or 
multidimensional representation, for example in three 
dimensions, is possible. 

[0028] In one particularly advantageous embodiment of the 
method, two locally coded magnetic resonance signals 
are combined such that the contribution of one magnetic 
resonance signal to the weighted combination is formed 
by multiplying this magnetic resonance signal by a 
weighting factor. The weighting factor depends on the 
other magnetic resonance signal such that it is greater 
for a large magnetic resonance signal than for a small 
magnetic resonance signal. 

[0029] In one specific embodiment, the weighting factor 
has a nonlinear dependency on the absolute value of a 
magnetic resonance signal . At least in one dimension of 
the raw magnetic resonance data, the same nonlinear 
dependency is used in order to have the same effect on 
the image sharpness, for example, over all the pixels 
in this dimension. 



New Patent Application 
Docket NO. 32860-000619/US 



[0030] In one particular embodiment of the method, the 
raw magnetic resonance data are processed using a 
plurality of filters before the differently filtered 
magnetic resonance signals are subsequently combined by 
means of weighted combination to form a single 
postprocessed magnetic resonance signal. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0031] A plurality of exemplary embodiments of the 
method, of the magnetic resonance tomography unit and 
of the magnetic resonance spectroscopy unit based on 
the invention will become more fully understood from 
the detailed description of preferred embodiments given 
hereinbelow and the accompanying drawings, which are 
given by way of illustration only and thus are not 
limitative of the present invention, and wherein: 

figure 1 shows a flowchart of an exemplary embodiment 
of the method which has been implemented in a 
magnetic resonance tomography unit, 

figure 2 shows a flowchart of an exemplary embodiment 
of the method in which two filtered magnetic 
resonance signals are subjected to weighted 
combination , 

figure 3 shows a numerical simulation of the 
processing of raw magnetic resonance data 
using a spin echo signal which is brought 
about by a square object in the image 
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gradient of a magnetic resonance tomography 
unit, 

figure 4 shows a numerical simulation of the 
processing when a low pass filter is applied 
to the spin echo signal from figure 3, 

figure 5 shows a numerical simulation of the 
processing when a high pass filter is applied 
to the spin echo signal from figure 3, 

figure 6 shows an abscissa detail from figures 3 to 5 
to illustrate the action of the method using 
the edge- emphasized profile which is achieved 
by the weighted combination, 

figure 7 shows an enlargement of figure 6 in the area 
of an edge of the square object, 

figure 8 shows a numerical simulation of the 
processing of raw magnetic resonance data for 
the case of magnetic resonance spectroscopy 
using an FID signal which has been formed by 
the sum of three Lorentz lines, and 

figure 9 shows a comparison of the magnetic resonance 
spectrum obtained by weighted combination 
with the low pass filtered magnetic resonance 
spectrum, both magnetic resonance spectra 
having been normalized to the same noise. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0032] Figure 1 schematically shows the flow of the 
method and the necessary components which are required 
within a magnetic resonance tomography unit 1 for the 
flow of the method. In a magnetic resonance tomography 
unit 1 having a conventional magnet and gradient system 
(not shown explicitly) used for local coding, a radio- 
frequency transmitter 3 is used to irradiate an 
examined object 5 with radio-frequency pulses. The 
electromagnetic signals emitted by the examined object 
5 are received using a radio- frequency receiver 7. The 
output signals from the radio- frequency receiver 7 form 
the raw magnetic resonance data in the form of a spin 
echo signal 9 . 

[0033] These data are postprocessed within a data 
processing installation 11. The latter subjects them to 
weighted combination with one another either in 
unfiltered form or in a form filtered using various 
filters 13 in a weighting unit 15. The result is a 
postprocessed magnetic resonance signal 17 in the form 
of a magnetic resonance picture. 

[0034] Figure 2 shows a flow chart of an exemplary 
embodiment of the method in which two filtered magnetic 
resonance signals are subjected to weighted combination 
in the data processing installation 11. In this case, 
the spin echo signal 9 is first demodulated in a 
demodulation unit 19, so that it is available in the k 
domain with the dimension of the sampling. Next, the 
spin echo signal 9 is supplied to two data processing 
units 21 having a respective filter. There, the spin 
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echo signal is filtered once using the low pass filter 
processor 23 and the other time using the high pass 
filter processor 25. 

[0035] Following the subsequent Fourier transformations 
27 and the absolute value formations 29, one firstly 
obtains the low pass filtered magnetic resonance 
picture 31 and the high pass filtered magnetic 
resonance picture 33, respectively. These two magnetic 
resonance pictures are combined in the weighting unit 
15 by weighted addition to form a single postprocessed 
magnetic resonance picture 35. The weighted addition is 
carried out by inputting parameters 37 using an input 
unit 39. The parameters 37 determine the weighting of 
the pixel 41 in the high pass filtered magnetic 
resonance picture 33 on the basis of the same pixel 43 
in the low pass filtered magnetic resonance picture 31. 
With a suitable choice of weighted addition and its 
parameters, the postprocessed magnetic resonance 
picture 35 has both a high SNR and great edge 
sharpness . 

[0036] When the method is used with a magnetic resonance 
spectroscopy unit, the flowchart is comparable with 
that in figures 1 and 2. The only difference is that 
raw data from magnetic resonance spectra are filtered 

and subjected to weighted addition by means of the 
method in order to achieve a higher resolution and a 
greater SNR. 

[0037] Figures 3 to 5 show the results of a numerical 
simulation of a spin echo signal for a square object in 
the field gradient of the magnetic resonance tomography 
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unit 1. The simulation relates to a section through a 
square object and illustrates the effect of the method 
on the object's edge sharpness. The object is in the 
area of the 128 pixels in the simulation. 

[0038] The figures respectively show the real and 
imaginary parts 51 and 53 of the echo signal after a 
filter has been applied to the echo signal. They show 
the filter function and the complex-value spin echo 
signal above a kx axis in the k domain. Various filter 
functions clarify the action of the filters on the echo 
signal in the K domain. In addition, figures 3 to 7 
show the result of Fourier transforming the filtered 
complex-value echo signal with subsequent absolute 
value formation. In this way, a simulated measurement 
result is obtained from a magnetic resonance tomography 
unit, in this case a one-dimensional section through 
the square object. The sectional view of the 
measurement result is plotted above the corresponding 
spatial axis. 

[0039] Figure 3 shows the real part 51 and the imaginary 
part 53 of the unfiltered echo signal and also the 
result of the Fourier transformation and of the 
absolute value formation for the unfiltered echo 
signal, in this case referred to as identity simulation 
result 57. In the simulation, no filter 55 was applied 
to the echo signal in this case. The oscillations in 
the identity simulation result 57 in the area of the 
rising and falling edges can be seen clearly. These 
oscillations limit the edge sharpness in the simulated 
(unpostprocessed) measurement. 
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[0040] Figure 4 shows the influence of the low pass raw 
data filter 59 on the spin echo signal. The low pass 
raw data filter 59 has been applied in the form of a 
Manning filter Fh to the spin echo signal, the Manning 
filter Fh having the following form in the 
one-dimensional k domain: 

Fj.^ =F„ =l/2yi[l + cos{2n{j -N/2)/N)] [Equation 2] 

[0041] In this case, N is the number of sampling points 
in the picture in one dimension. The filter modulation 
follows a cosine function whose maximum is shifted by 
N/2, that is to say to the center of the k domain. The 
index j describes the jth sampling point- The filter 
suppresses the marginal area of the k domain with lower 
intensity in the echo signal. This can clearly be seen 
from the reduced oscillation amplitude of the filtered 
spin echo signal in the marginal area. The low pass 
simulation result 61 in figure 4 shows a rounded edge 
for the square object. This can also be seen in figures 
5 and 6 . 

[0042] In figure 5, the simulation is carried out for a 
high pass raw data filter 63, with the high pass filter 
in the k domain being a modulated Hamming filter Fh 
provided with an offset: 

Fh;h = (Fh + Offset)^, cos" {tt {j - N/ 2)/N) [Equation 3 ] 

[0043] In this case, the modulation of the Manning filter 
Fh takes place with a cosine function to the power of n, 
with the minima of the cosine function being at the two 
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ends and in the center of the k domain. The additional 
minima at the ends have the advantage that they cause 
no additional artefacts in the Fourier transformation. 
Thus, by way of example, the "truncation artefact" in 
the Fourier transformation is suppressed. 

[0044] In this case, the offset parameter of the 
modulated Hanning filter has the value 0.2 and the 
power n of the cosine modulation has the value 1 . The 
high pass raw data filter 63 amplifies the wings of the 
echo signal in the k domain. In the real part 51 and in 
the imaginary part 53 of the filtered echo signal, it 
is possible to see the amplified amplitude of the echo 
signal oscillations in the marginal areas and a 
suppression in the central area. The high pass 
simulation result 65 shows a marked overshoot in the 
area of the edges of the square object. 

[0045] Figure 6 shows a detail from the simulation 
results in the area of the pixels 35 to 55, i.e. in the 
area of the rising edge. The oscillations in the 
identity simulation result 57 can clearly be seen. The 
edge of the square object is rounded in the case of the 
low pass simulation result 61. 

[0046] The high pass simulation result 65 shows a 
pronounced signal overshoot right in the area of the 
rise. Using the weighted addition of the two filtered 
simulation results, it is now possible to achieve an 
edge - emphasized postprocessing result 67. To this 
end, a weighted combination has been performed in line 
with the equation 1 already indicated above. In this 
case, the low pass simulation result 61 corresponds to 
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the magnetic resonance signal A and the high pass 
simulation result 65 corresponds to the magnetic 
resonance signal B. 

[0047] The parameters X and K are set such that a maximum 
sharpness is produced in the edge area. In this case, X 
and K are usually between the values 1 and 3. The level 
obtained for the weighting of the contribution of B 
depends on the ratio of the absolute value in the 
respective pixel in B to the maximum value Anvax of the 
pixels in A, i.e. the higher the value of the high pass 
simulation result 65 in a pixel, the higher its 
contribution. X then performs a type of linear 
weighting, while K weights the influence of this 
quotient as a power of the ratio of B and A, that is to 
say indicates how greatly the correction of A is 
determined by the ratio of magnitudes for B with 
respect to A. The parameters X and K have been set to 
the values 0.6 and 2.3, respectively, in figures 6 and 
7 . 

[0048] Besides the edge sharpness, the SNR of the 
postprocessing result 67 is also increased, and is now 
comparable to the SNR of the low pass simulation result 
61. This can be seen from the suppression of the signal 
before the rising edge, i.e. there is a low base 
contribution outside the object. 

[0049] Figure 7 shows a further enlarged detail from the 
simulated measurement results . The representation 
focuses on the edge area of the square object in the 
area of pixels 42 to 55 with signal levels between 10 
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and 13. It can clearly be seen how the greatly 
oscillating identity simulation result 57 is freed from 
the radio- frequency oscillations using the low pass 
filter. Next, the weighted addition fills the corner 
area of the square object on account of the high signal 
absolute values in the high pass simulation result 65. 
The oscillations in the identity simulation result 57 
are also called "ringing". This artefact, which is also 
called truncation artefact, is suppressed using the 
method of the invention. 

[0050] Figure 8 shows the flow of the method applied to 
raw magnetic resonance data from a simulated magnetic 
resonance spectroscopy measurement in which an FID 
signal 71 is measured which is formed by the sum of 
three Lorentz lines (73, 75, 77). Figure 8 shows the 
real parts of the FFT spectra involved in the method, 
in the simulation, a noise contribution has been added 
to the synthesized FID signal 71. Low pass filtering of 
the FID signal 71 results in the low pass filtered 
spectrum 79, and high pass filtering results in the 
high pass filtered spectrum 81. Following the weighted 
combination of the two filtered spectra, the 
postprocessed magnetic resonance spectrum 83 is 
obtained. It is possible to see markedly reduced noise 
and high resolution of the three Lorentz lines on 
account of the use of the method of the invention. 

[0051] In figure 9, the postprocessed magnetic resonance 
spectr^am 83 is compared with the low pass filtered 
spectriam 79. To clarify the influence of the method, 
the two spectra have been normalized to the same noise 
in this case. The narrower line widths and a 
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considerable increase in the contrast in the 
postprocessed magnetic resonance spectrum 83 are 
clearly evident . 

[0052] Exemplary embodiments being thus described, it 
will be obvious that the same may be varied in many 
ways. Such variations are not to be regarded as a 
departure from the spirit and scope of the present 
invention, and all such modifications as would be 
obvious to one skilled in the art are intended to be 
included within the scope of the following claims. 



19 



